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Abstract 

The higher twist corrections h N (x)/Q 2 to the spin dependent proton and neutron 
<7i structure functions are extracted from the world data on gi(x, Q 2 ) in a model 
independent way and found to be non-negligible. Their role in determining the 
polarized parton densities in the nucleon is discussed. It is also considered how 
the results are influenced by the recent JLab and HERMES / d inclusive DIS data. 
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1 Introduction 



Spurred on by the famous EMC experiment [1] at CERN in 1987, there has been a huge 
growth of interest in polarized DIS experiments which yield more refined information 
about the partonic structure of the nucleon, i.e., how the nucleon spin is divided up 
among its constituents, quarks and gluons. Many experiments have been carried out 
at SLAC, CERN, DESY and JLab to measure the longitudinal (A\\) and transverse 
(A±) asymmetries and to extract from them the photon- nucleon asymmetries Ai(x, Q 2 ) 
and A 2 (x, Q 2 ) as well as the nucleon spin-dependent structure functions g±(x, Q 2 ) and 
g 2 (x,Q 2 ). 

There is, however, an important difference between the kinematic regions of the 
unpolarized and polarized data sets. While in the unpolarized case we can cut the low 
Q 2 and W 2 data in order to eliminate the less known non-perturbative higher twist 
effects, it is impossible to perform such a procedure for the present data on the spin- 
dependent structure functions without loosing too much information. So, to extract 
correctly the polarized parton densities from the experimental data a special attention 
should be paid to the higher twist (powers in l/Q 2 ) corrections to the nucleon structure 
functions. Their role in determining the polarized parton densities in the nucleon 
using different approaches of QCD fits to the data is discussed in this talk. It is also 
considered how the results are influenced by the recent JLab [2] and HERMES/d [3] 
data. 

2 QCD treatment of gi(x,Q 2 ) 

In QCD the spin structure function gi can be written in the following form: 

g x (x, Q 2 ) = gi(x, Q 2 ) LT + g x (x, Q 2 )m , (1) 

where "LT" denotes the leading twist (r = 2) contribution to gi, while "HT" denotes 
the contribution to g\ arising from QCD operators of higher twist, namely r > 3. In 
(1) we have dropped the nucleon target label N. The HT power corrections (up to 
0(1/ Q 2 ) terms) can be divided in two parts: 

gi(x, Q 2 )ht = h(x, Q 2 )/Q 2 + h TMC (x, Q 2 )/Q 2 , (2) 

where h™ c (x, Q 2 ) are the exactly calculable [4] kinematic target mass corrections and 
h(x, Q 2 ) are the dynamical higher twist (r = 3 and r = 4) corrections to g±, which are 
related to multi-parton correlations in the nucleon. The latter are non-perturbative 
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effects and cannot be calculated without using models. gi(x, Q 2 )lt in (1) is the well 
known pQCD expression and in NLO has the form 

9i(x, Q 2 ) pQ cd = \ £ eli(A q + Aq) ® (1 + ^f±5C q ) + ® ^], (3) 

where Aq(x,Q 2 ), Aq(x,Q 2 ) and AG(x,Q 2 ) are quark, anti-quark and gluon polarized 
densities in the proton, which evolve in Q 2 according to the spin-dependent NLO 
DGLAP equations. 5C(x) q: o are the NLO spin-dependent Wilson coefficient functions 
and the symbol <S> denotes the usual convolution in Bjorken x space. Nf is the number 
of active flavors. 



3 QCD fits to the data and the role of higher twists 

Up to now, two approaches have been mainly used to extract the polarized parton 
densities (PPD) from the world polarized DIS data. According to the first [5, 6] the 
leading twist LO/NLO QCD expressions for the structure functions g± and F^ have 
been used in order to confront the data on Ai(fa gi/Fi) and g\jF\. It was shown [6, 7] 
that in this case the extracted from the world data "effective" HT corrections h Al (x) 
to A 1 

are negligible and consistent with zero within the errors, h Al (x) (see Fig.l). This 
result has been confirmed independently in [5]. In Fig. 1 are also shown (open circles) 
our new results on the HT corrections to A 1 including in the world data set [1, 8] the 
recent JLab [2] and HERMES [3] data on g\jF x for neutron and deuteron, respectively. 
As seen from Fig. 1, due to the much more precise JLab and HERMES new data, the 
HT corrections h Al (x) to Ai for the neutron and deuteron targets are much better 
determined now at large x and better consistent with zero in this kinematic region. 

What follows from these results is that the higher twist corrections to g\ and F\ 
compensate each other in the ratio g\jF\ and the PPD extracted this way are less 
sensitive to higher twist effects. 

According to the second approach [9], g\jF x and A x data have been fitted using 
phenomenological parametrizations of the experimental data for F 2 (x, Q 2 ) and R(x, Q 2 ) 
(Fi has been replaced by the usually extracted from unpolarized DIS experiments F 2 
and R). Note that such a procedure is equivalent to a fit to (gi) eX p, but it is more precise 
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than the fit to the g\ data themselves actually presented by the experimental groups. 
The point is that most of the experimental data on g\ have been extracted from the A\ 
and gi/Fi data using the additional assumption that the ratio gi/Fi does not depend 
on Q 2 . Also, different experimental groups have used different parametrizations for F 2 
and R. 

If the second approach is applied to the data, the "effective higher twist" contribu- 
tion h Al (x)/Q 2 to Ai(gi/Fi) is found [5] to be sizeable and important in the fit [the 
HT corrections to g\ cannot be compensate because the HT corrections to F 1 (F 2 and 
R) are absorbed by the phenomenological parametrizations of the data on F 2 and R\. 
Therefore, to extract correctly the polarized parton densities from the g 1 data, the HT 
corrections to g x have to be taken into account. Note that a QCD fit to the data in this 
case, keeping in gi(x, Q 2 )qod only the leading-twist expression (as it was done in [9]), 
leads to some "effective" parton densities which involve in themselves the HT effects 
and therefore, are not quite correct. 

Keeping in mind the discussion above we have analyzed the world data [1, 8] on 
inclusive polarized DIS taking into account the higher twist corrections to the nu- 
cleon structure function g^(x : Q 2 ). In our fit to the data we have used the following 
expressions for gi/F 1 and A^. 

^f(x,g 2 )] g?(x,Q 2 ) LT + h N (x)/Q 2 n [l + R(x,Q 2 ) exp ] 

„r,N zx- 



F 2 N (x,Q 2 ) exp ~ (1 + 7 2 ) 



Ffix^Q 2 ) 

A»(x,Q 2 ) exp # 9 ' {X > g^ LT + 2 f {X)/Q2 2x[l + R(x, Q 2 ) exp ] , (5) 

r 2 [X, j exp 

where gi(x, Q 2 )lt is given by the leading twist expression (3) including the target mass 
corrections (N=p, n, d). The dynamical HT corrections h N (x) in (5) are included and 
extracted in a model independent way. In our analysis their Q 2 dependence is ne- 
glected. It is small and the accuracy of the present data does not allow to determine 
it. For the unpolarized structure functions F^(x, Q 2 ) exp and R(x, Q 2 ) eX p we have used 
the NMC parametrization [10] and the SLAC parametrization R4998 [11], respectively. 
The details of our analysis are given in [12]. We have found that the fit to the data 
is significantly improved when the higher twist corrections to g\ are included in the 
analysis, especially in the LO QCD case. We have also found that the size of the HT 
corrections to g\ is not negligible and their shape depends on the target (see Fig. 2). 
In Fig. 2 are also presented (open circles) our new results on the HT corrections to 
gi including in the world data set the recent JLab [2] and HERMES [3] data. As seen 
from Fig. 2, the higher twist corrections to the neutron spin structure functions in the 
large x region are much better determined now. It was also shown (see Fig. 3) that 
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the NLO QCD polarized PD(^ T + HT) determined from the data on g>i, including 
higher twist effects, are in good agreement with the polarized PD(gf LO / F± LO ) found 
earlier from our analysis [6] of the data on gi/F\ and A\ using for the structure func- 
tions g± and F\ only their leading twist expressions in NLO QCD. This observation 
confirms once more that the higher twist corrections to g±/F 1 and A 1 are negligible, 
so that in the analysis of gi/F 1 and A 1 data it is enough to account only for the lead- 
ing twist of the structure functions g± and F\. On the other hand, in fits to the g\ 
data themselves the higher twist contribution to g\ must be taken into account. The 
latter is especially important for the LO QCD analysis of the inclusive and SIDIS data. 

This research was supported by the JINR-Bulgaria Collaborative Grant, by the 
RFBR (No 02-01-00601, 03-02-16816), INTAS 2000 (No 587) and by the Bulgarian 
National Science Foundation under Contract Ph-1010. 
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Figure Captions 



Fig.l. Effective higher twist contribution h Al (x) to the spin asymmetry A±(x, Q 2 ) 
extracted from the data. 

Fig. 2. Higher twist corrections to the proton and neutron g 1 structure functions 
extracted from the data on g 1 in the NLO QCD approximation for gi(x, Q 2 )lt- 

Fig. 3. NLO(JET) polarized parton densities PD(#f LO + HT) (solid curves) together 
with error bands compared to PD(^ LO /F 1 NLO ) (dashed curves) at Q 2 = 1 GeV 2 . The 
error bands represent the total errors. 
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